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Abstract. Methane brown dwarfs {CHiBDs), i.e. ob- 
jects similar to the low mass star Gliese 229B, are char- 
acterized by very unusual colours: 3-K^O and I-J^5. An 
analysis of the ESO public images, which cover an area of 
37 sq. arcmin in the three filters, yields one clear CH^BD 
candidate with J=20.2 and two fainter sources with J~23. 
The resulting observed density of methane dwarfs ranges 
from R:!lOO to «500 objects per sq. degree with J<23.1, 
quite high but compatible, within the errors, with what 
we derive from simulations of the stellar population of 
brown dwarfs in the disk of the Galaxy adopting reason- 
able (although still speculative) hypotheses on the evolu- 
tion of BD colours with cooling. The predictions presented 
here can be useful in constraining the results from future 
searches of these objects in the infrared. Deep imaging of 
several sq. degrees of sky down to J~22, K~22 and 1^27 
are necessary to significantly improve the observational 
database. 



1. Introduction 

Until recently, only three very low luminosity objects were 
known with spectral type later than M9.5, but the two 
near IR ground based surveys DENIS and 2MASS opened 
a new window in the study of the end of the main sequence 
and its continuation into the realm of Brown Dwarfs 
(BDs). In fact, DENIS added three objects (Delfosse et 
al. |1999| Tinney et al. |1998D , and 2MASS (Kirkpatrick et 
al. 1998) increased this sample till now by about 20 ob- 



jects spectroscopically confirmed by the Low Resolution 
Infrared Spectrometer at Keck. These findings are defining 
a new spectral class, the "L dwarfs" spanning Toff~2200- 
1500 K, their exact temperatures depending on the treat- 
ment of dust in model atmospheres. 

It was recently recognized that the near IR colours of 
the late M dwarfs (Toff^2600 K) could not be reproduced 
by model atmospheres, unless they included consideration 
of grains. Although the spectral feature are very com- 
plex, the first generation of model atmospheres including 



grain condensations (Tsuji et al. 1996, AUard 1998] ) are 
able to interpret the main emerging spectral characteris- 
tics, and include predictions also for the L-dwarf spectra. 
Dust grains forming in the upper layers of very cool dwarfs 
have a large opacity, and produce a redistribution of the 
IR flux (greenhouse effect), explaining the very red J-K 
colours of the latest M and of the L-dwarfs. The L-dwarfs 
clearly define the low end of the MS, which, depending on 
the opacity, and thus on the dust inclusion in the mod- 
els, m ay be shifted even as low as '^1800 K (Burrows and 
Sharp 1999). However, due to the steepness of the mass 



luminosity relation around the hydrogen burning mini- 
mum mass (DAntona & Mazzitelli 1985), the range of 
Tcff's of L-dwarfs should be mostly populated by cooling 
brown dwarfs. An indication confirming this interpreta- 
tion is the presence of Lithium in about half of the 2MASS 
objects, indicating that they are unable to burn their ini- 
tial lithium, and thus that their mass is below ~0.06 Mq 



(see however Allard 1998 for a different interpretation). 

When L-dwarfs evolve further, they will eventually 
reach Teff<1500 K. Just as for normal white dwarfs, the 
evolutionary time is dictated by the release of the thermal 
energy of the ions, so that it increases when decreasing 
the luminosity. Only one object was until recently iden- 
tified as the result of such cooling: the companion to the 
low mass star Gliese 229, having a spectrum showing the 
CH4 bands, predicted by theory. This object, the "Rosetta 
stone" of cool BDs, at Toff ~900 K and luminosity well be- 
low 1O~^L0 (Leggett et al. 1999), shows in fact a peculiar 
spectrum characterized by metal depletions, and domi- 
nated by H2O, CII4 and alkali metals. The ffux redistri- 
bution at very low T^s produces a spectrum much bluer in 
J"K, while keeping a very red I-K colour (e.g. Burrows et 
al. 1997). Therefore, while the typical colour of L-dwarfs 
are J-K>1.3 (e.g. Leggett et al. 1998), the cooler BDs 
should have J-K^O, as confirmed by the colour of Gliese 
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229B (Nakajima et al. |1994 Legget et al. |1999 ). 

Although this outline of low mass and brown dwarf 
evolution is qualitatively sound, we are quite far from a 
complete theoretical understanding of the problem. At the 
upper mass end (M ^ 0.4Mq) the inclusion of molecular 
opacities in the model atmospheres, and the use of bound- 
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Fig. 1. HR diagram with the evolutionary tracks of low 
mass stars. The dashed lines are from D'Antona & Mazz- 



itelli ( |1997|) and cover the 0.6 to 0.2 Mq range, while the 
tracks for lower mass stars (0.237 to 0.03 Mq) are from 
Burrows et al. (199J). The square symbols mark the po- 
sitions at 10^" years. The approximate boundaries of the 
M, L and methane dwarfs are shown. The Main Sequence 
limit is between 2000 and 1800K depending on the opac- 
ities adopted, so it falls in the middle of the L spectral 
type. 



ary conditions from non-gray model atmosphere for the 
interior models has improved considerably the match be- 
tween the stellar models and the observed colour of stars 



(see Baraffe et al. 1998, hereafter B98 ). We have remarked 
that the optical colours at Teff<2500 K require grain in- 
clusion in the models. This is a difficult task, although 
insig ht is e xpected based on the "dusty" modeling by Al- 
lard (1998k and on the grain clouds approach by Burrows 
et al. ( 1998| ). Nevertheless, the unknowns in the modeling 
of grains (e.g. the particle size distribution, the spatial dis- 
tribution of clouds, their composition, or the assumption 
of chemical equilibrium without rain-out) shift very far 
away in time a self-consistent approach. In addition, the 
detailed colour transition from L dwarfs to CH^EDs can 
not be yet described quantitatively. 

Fig. 1^ shows the evolutionary tracks of dwarfs and 
BDs, and a very schematic separation between M, L and 
methane dwarfs, in terms of Tcff , which we choose accord- 
ing to Reid et al. 1998. These limits must be taken with 



the caution due to the uncertainties described above. 
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Fig. 2. Selected regions of the ESO deep fields showing 
the CH4BDS candidates. Each sub- frame is 28"x28" with 
N to the left and E to the bottom, the position of the can- 
didate is marked by a white rectangle. The object magni- 
tudes and colours are ±0.15 while the upper limits are at 
the 3cr level. 



2. Observational constraints on the density of 
faint CH4.BDS 

In spite of all the difhculties in the models, the very blue 
J-K colour, together with a very red I-K, should be a 
signature of the very low Teff, older BDs which we will 
call CH4BDS. This simple recognition has prompted us 
to search for faint objects characterized by small J-K but 
very large I-K or R-J, indicative of low Toff. A very useful 
data set comes from the ESO public survey data which 
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Table 1. Adopted colours and bolometric corrections 



(1) 



(2) 



(3) 



logTeft BCy 



(1) 



(J-K)!^L 






2.768 


2.23" 


-1.44" 


-1.44 


3.000 


2.20" 


-0.10" 


-0.10 


3.301 


2.19 


0.87 


1.37 


3.362 


2.16 


0.92 


1.21 


3.423 


2.07 


0.93 


1.06 


3.447 


2.00 


0.92 


0.97 


3.498 


1.88 


0.87 


0.91 


3.529 


1.81 


0.84 


0.84 


3.554 


1.75 


0.81 


0.81 


3.592 


1.66 


0.79 


0.79 


3.631 


1.54 


0.71 


0.71 


3.674 


1.39 


0.57 


0.57 


3.713 


1.26 


0.43 


0.43 



Bol ometr ic correction in J, from the models of Burrows et 
al. (1998) except for those in note a 



Colours predicted by models, from B98 except for those in 

note a 

Empirical colours adopted h ere, s ee Sect, p.^ 

Values from Burrows et al. (1997) 
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Fig. 3. Example of the predicted distribution of apparent 
stellar J magnitudes versus the distance, for simulated ob- 
servations covering 1 square degree. 



are targeted at regions far from the galactic plane and in- 
clude the N TT SUSI deep fiel d (^= 284°, 6=+53°, Arnouts 



199c, Moorwood et al. 1998), the EIS HDF-S field 



-49°, da Costa et al. |1999D and the A XAF 
=224° 



et al. 

(^=328°, 

deep field (^=224°, 6=-54°, Rengelink et al 
overlap between infrared (K, J) and optical (I, R, V) im- 
ages cover a total area of 37 sq. arcmin and is complete to 
K=22.6, J=24.4 and 1=26.3 mag. These limits are some- 
what unbalanced for our purposes, i.e. too shallow in K 
and I for the depth reached in J. The colour selection we 
use to extract methane dwarf candidates effectively trans- 
lates into a cut at J=23.1 and much above the limit of 
completeness of the J images. 

The methane dwarf candidates were selected as those 
point-like objects with J-K<0.5 and I-J>3, i.e. colours 
which are incompatible with any known stellar source 
other than GL229B. This yielded three objects whose find- 
ing charts and magnitudes are shown in Fig. g. Considered 
the statistical errors, the corresponding density ranges be- 
tween «100 and «500 objects per sq. degree with magni- 
tudes J<23.1. 

In should be noted that the colours of the selected 
objects are also compatible with quasars or QSOs at 
Z'^9, i.e. at redshifts where the Lya falls into the J band 
and the Lyman-break is moved beyond the I photometric 
band. In such a case an observed magnitude J=23 would 
translate into an absolute blue magnitude Mb — —25 
and typical of low-medium luminosity QSOs. However, 
we consider this possibility unlikely because current 
estimates of the density of quasars at various redshifts 



indicate a rapid drop beyond z—3 with basically no 
object at z > 6 (e.g. Shaver et al. 1996] ). Encouragingly, 
the CH4BD nature of the brightest object was recently 
confirmed by spectroscopic IR observations at the VLT 



19991) . The (Cuby et al. |I999|) 



3. Modelling of number counts and predicted 
colours 

As outlined in the introduction, it is not possible to model 
self-consistently the whole BD evolution. We decided 
therefore to adopt an admittedly rather simple-minded 
approach to the colour derivation, in order get a first es- 
timate of the expected number counts of CH4BDS. On 
the other hand, we took care to model carefully the galac- 
tic disk, so that the global model can be easily improved 
when more precise evolutionary tracks will be available. 

3.1. The theoretical models 

As a first input, we need the time evolution of logL/LQ 
and Toff for stellar masses which trace the low end of the 
main sequence and the brown dwarf evolution. We will 



use as a basis the models by Burrows et al. (1997 and 



private communication) employing non gray model atmo- 
spheres up to 1300K, and similarly low resolution non 
gray models with grains up to 2400K. Above this temper- 
ature, gray models with Rosseland opacities from Saumon 
are adopted. These tracks are computed for masses from 
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Fig. 4. Stellar counts (in bins of 0.5mag, per square de- 
gree) as function of the apparent J magnitude for three 
selected mass ranges. 



0.237Mq to O.O3M0, and are evolved up to an age of 
20Gyr, which is reached at quite low T^s for the smaller 
masses, so they constitute a very complete set to deal with 
all the very low mass and brown dwarf regimes, from the 
M spectral type, through the L dwarfs, to the postulated 
CH4BDS phase. Above this mass, for the purpose of hav- 
ing predictions available also for somewhat hotter stars, 
we adopt the tracks by D'Antona and Mazzitelli ( 1997 ) 
from 0.25 to 0.8Mq. Although the latter are computed 
with somewhat different inputs (e.g the deuterium initial 
abundance is a half of the abundance assumed by Burrows 
et al., the treatment of over-adiabatic convection is differ- 
ent) the merging of the two sets is very reasonable, also 
because mainly the main sequence or cooling part of the 
sequences are important to determine the stellar counts. 
The gray atmosphere low mass models by D'Antona and 
Mazzitelli (1997) predict luminosities larger by ~ 15% 
(generally <0.1mag) than the non gray models e.g. by 



Chabrier & Baraffe (1997, hereafter pB97| ) and have T^e 
~ 2% larger. The Burrows et al models above 2400K com- 
pare well with those of |CB97] in luminosity, and they differ 
by 2% in Teff. 

3.2. Colours and bolometric corrections 

The magnitudes in I, J and K are not directly available in 
the track sets, so we decided to adopt the admittedly crude 
procedure of estimating the I and J band bolometric cor- 
rections (BC) from the B98| mo dels, which are consistent 
with the |CB97 models. The B9^ values have been used for 



Fig. 5. Stellar counts in bins of half magnitude and per 
square degree as a function of the apparent magnitude J, 
in the case 1=0, b=90, x=0. The cumulative counts from 
0.03 to O.8M0are shown in the thick line. Partial counts 
due to separate mass ranges are also given to show at 
which magnitudes each mass range contributes mostly to 
the counts. The continuous line at the lower right repre- 
sents the counts from 0.03 to O.O4Af0, the dash-dotted 
one from 0.03 to O.O8Af0, the dotted line from 0.03 to 
O.12M0and the dashed from 0.03 to 0.2MoThe BD mass 
range is largely represented only for J>20. The counts are 
limited to J=25, but BDs would be present also above this 
magnitude. 

rcff> 2000 K. We extrapolate these BCs to Teff=600 K 



using the models of Burrows et al. (1997). We remark that 
the J band counts should be little affected by uncertain- 
ties in the bolometric corrections. In fact, BCj is 2.19mag 
at 2000K, and very similar (2.23mag) at Teff=600 K (see 
table |l|) . In other words, 25-30% of the flux for these ob- 
jects emerges in the J band, regardless of their tempera- 



ture (e.g. Burrows et al. 1998) 



The colours J-K are taken from Burrows et al. 1998 for 



rcff<1000 K. Above 3300 K we adopt the colours of B98 



In the intermediate range, these models fail reproduce the 
observations of cooler dwarfs (Leggett et al. 1998[ |B98[ ) 
and computations including grains should be used (e.g. 
Allard |1998| Jones & Tsuji |1998| ) to get the very red J-K 
colours down to GD165B. However, dusty models are still 
very incomplete and we therefore employed a semiempir- 
ical approximation, imposing the models to pass through 
the colour of GD165B and then to evolve to the blue. 
Table |l| shows the adopted relation between Tcs and J- 
K. The reader should be well aware that this approach 
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represents an educated guess on the CH4BDS colour evo- 
lution. Nevertheless the semiempirical approach mainly 
influences the relative number counts of the J versus K 
band, but not the Tofi selection, which however is itself is 
uncertain in the absence of better modeling. 

3.3. The integrated number counts 

The number of stars within a given range of temperature 
(TeffzbAT) and apparent magnitude (.T +A .T) expected from 
observations directed at a given direction I, b (galactic co- 
ordinates) and covering a relatively small solid angle Q is 
given by 



10 



/•OO 

E(Teff,J)= / $[Teff,(J-51ogi?io)]p(i?)r!i?2di? 
Jo 



(1) 



where R is the distance from the sun along the I, b line 
of sight, i?io = (-R/10 pc) and p is the total stellar density 
whose radial dependence can be approximated by (Bachall 
& Soneira, (igTSl) 



p{R) 



Po exp{-z/H - (r - ro)/h) 

R sin 6 

■\/ {R cos 6)2 — 2Rro cos b cos I 



The coordinate z is measured perpendicular to the galactic 
plane, r is the radial distance from the galactic center and 
ro=8 kpc is the distance of the Sun. The scale height is 
if =325 pc, consistent with the counts of M dwarfs, and 
the scale length is h^3.5 kpc. The local stellar density po 
is given by 



Po = 0.336 



AU 



M„ 



M 



0.1 M, 



^(l+x) 



dM 



where we have adopted a power law initial mass func- 
tion dN/dM oc M~'^^+^^ and the norma lization factor at 
0.1 Mq is taken from Reid & Gizis ( 1997 ). The mass limits 
considered here are Mmm=0.03 Mq and Mmax=Q-S> Mq. 
In Eq. (1) the term $(Teff,Mj) gives the fraction of 
stars within the selected range of temperature and ab- 
solute J magnitude, i.e. the relative density of stars on 
a given position of the theoretical Teff-Mj HR diagram. 
This function, which depends on the shape of the IMF and 
on star formation history, is constructed via a Montecarlo 
technique assuming an uniform birth-rate from t—0 to the 
present disk age idisfe = 10 Gyr. The program randomly ex- 
tracts masses (with probability given by the chosen IMF 
shape) and ages from the track sets of Fig. |l|. Each ex- 
traction provides a point in the Mboi~Tcff diagram which 
translates into absolute magnitudes using the bolometric 
corrections of Table El. The extracted stars are grouped 
into suitably small bins of temperatures and Mj, this 
yields the matrix which represents the function $ which 
is used to compute the integral of Eq. (1). 
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Fig. 6. Example of the predicted distribution of apparent 
stellar J magnitudes versus J-K colour for observations 
covering of 1 square degree. The points cluster along the 
main sequence defined by the thick part along J-K~l, 
then we recognize the L-dwarfs which populate the re- 
gion at 1<J-K<1.6, before evolving to much bluer colours 
when the methane bands begin to appear in the atmo- 
sphere. 



3.4. Results 

The distribution of the apparent magnitude J as a function 
of the distance (Fig. 0) shows that the bulk of M dwarfs 
is at magnitudes J<20, where BDs begin to appear. L- 
dwarfs are numerous at J^20-22, while CH4BDS begin 
to be common only at J>21. 

Fig. ^ and Fig. |^ also show the same feature. The first 
displays the number counts per J magnitude within three 
mass bins: dwarfs (0.5-0.7 Mq), low mass dwarfs (0.08- 
0.14 Mq) and brown dwarfs (below 0.08 Mq). Being rel- 
atively bright, the dwarfs show a maximum in the dis- 
tribution at J~17, the lowest mass dwarfs also show a 
maximum at J~20, while the BDs number increases down 
to the limit at which we count them (J=25). Therefore, 
either very deep searches are necessary to show an ap- 
preciable number of CH4BDS, or a relatively deep survey 
covering most of the sky. Fig. shows the total number 
counts versus J magnitude, putting into evidence which 
mass ranges contribute to its different parts. The broad 
maximum in the counts at J=17-19 includes only dwarfs 
and results from the combination of the effects due to the 
IMF, the mass luminosity relation for the main sequence 
and the radial density distribution. The rising of the num- 
ber counts at J>22 is instead due to the low mass brown 



F. D'Antona E. Oliva, and A. Zeppieri: The density of CH4 Brown Dwarfs 



Table 2. Predicted brown dwarf densities*^^-' 
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Z=0° 6=90 
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14.0 


.0004 + .0003 


.004 


.01 


.3 


7 


.0004 + .0003 


.004 


.01 


.3 


7 


15.0 


.001 + .001 


.016 


.05 


1 


20 


.001 + .001 


.016 


.05 


1 


25 


16.0 


.006 + .005 


.06 


.19 


4 


65 


.006 + .005 


.06 


.19 


5 


78 


17.0 


.022 + .019 


.20 


.65 


12 


160 


.023 + .020 


.21 


.80 


13 


215 


18.0 


.085 + .074 


.78 


2.1 


33 


320 


.088 + .079 


.84 


2.4 


40 


490 


19.0 


.32 + .27 


2.6 


6.4 


79 


515 


.33 + .30 


2.9 


7.6 


106 


900 


20.0 


1.2 + .91 


7.7 


16 


150 


650 


1.3 + 1.0 


9.1 


22 


240 


1300 


21.0 


4.0 + 2.7 


19 


34 


230 


702 


4.4 + 3.2 


25 


51 


420 


1500 


22.0 


13 + 7.5 


40 


54 


280 


710 


15 + 9.3 


60 


95 


570 


1550 


23.0 


38 + 16 


74 


73 


290 


710 


50 + 20 


115 


130 


570 


1550 


24.0 


90 + 30 


90 


74 


290 


710 


112 + 40 
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575 


1550 
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160 + 38 
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270 + 72 
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.0008 + .0006 


.006 


.018 


.3 
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.006 


.018 


.3 
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15.0 
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.025 


.068 


1 


9 
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.025 


.070 


1 


9 


16.0 


.01 + .009 


.1 


.25 


4 


32 
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.1 


.25 


4 


32 


17.0 


.045 + .036 


.31 


.9 


12 


68 
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.32 
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13 
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.17+ .13 


1.2 


3 


32 
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3 


40 
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240 
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4.5 


10 
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2.3 + 1.6 


12 


22 
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14 


28 
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600 


21.0 


8.0 + 4.9 


30 


44 


225 


345 


9.0 + 5.7 


39 


66 


400 


720 


22.0 


25 + 13 


61 


69 


270 


350 


29 + 17 


94 


120 


550 


760 


23.0 


84 + 24 


99 


73 
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350 
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550 


760 


24.0 


200 + 38 
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76 
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240 + 95 


235 
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560 


760 


25.0 


385 + 52 


130 


76 
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350 


535 + 130 
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(1) 



Cumulative number of objects per square degree with magnitude J<Jii 

reff<1400 K + 1400 K<Tcff<1600 K 

1600 K<Teff<2200 K 

2200 K<Teff<2500 K 

2500 K<Tcff<3000 K 

3000 K<Tcft<3500 K 

IMF in the form dN/dM oc M-(i+^) 



dwarfs. Finally, Fig. g shows simulated J vs J-K colour- 
magnitude diagram for observations covering 1 sq. degree 
of sky. The large scatter of J-K colours for the CH4BDS 
reflects the effect of methane bands absorption on the J- 
K colour, which is expected to increase rapidly between 
Toff=1000 K (J-K=-0.1) and reff=600 K (J-K=-1.4, see 
table |I|). 

The resulting counts are tabulated in Table g for galac- 
tic coordinates similar to those of the ESO fields, namely 
6=55° I =340°, and at the galactic pole (6=90°). 

3.5. Comparison with observations 

Table shows that a considerable number of L-dwarfs are 
predicted, but practically no CH4BDS (at most 1 in 1000 
sq. degrees), in a large area limited at K<14.5, in agree- 
ment with the early results of 2MASS. The very recent 
obser vation s of four relativ ely b right CH^BDs (Strauss 
et al. 1998| , Burgasser et al. 1999 ) yields a surface density 



of ~0.002 cool methane dwarfs with J<16 and a factor of 
3-6 lower than those predicted here using flat and steep 
IMFs, respectively, and considering only the density of 
objects of Tcfr<1400 K. However, these differences are not 
necessarily significant, given the small number of objects 
detected and the possibility that the the above searches 
may not be yet complete. 

Our observed density of 3 objects with J<23.1 within 
37 sq. arcmin is consistent (i.e. within the Poisson's noise) 
with the space density predicted using a relatively steep 
IMF {x — 1) while it is a factor of 4±2 larger than ex- 
pected from a flatter IMF {x = 0). Remember that this 
latter IMF is certainly more consistent with the slope de- 
rived from nearby low mass stars (Reid and Gizis 1997). 
Given the very small number of objects detected we con- 
sider it premature to draw strong conclusions about the 
IMF. 
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4. Conclusions 

We presented three CH^BDs candidates with J<23.1 
found in an area of 37 sq. arcmin covered by the ESO 
pubhc images survey and modelled the expected number 
counts for these objects as well as for L and M dwarfs. 
Theoretical predictions are made on the basis of a sound 
galactic model including the most complete set of today's 
available evolutionary tracks and colours-Toff correlations. 
The model could be significantly improved only when self- 
consistent computations will be available, including model 
atmospheres tested on the methane dwarfs which will be 
discovered. 

The predicted surface densities are in reasonably good 
agreement with the observed counts of L dwarfs and 
with the sparse observational data presently available for 
methane dwarfs, namely the lack of objects with K<14.5 
in the early 2MASS survey, the recent finding of four rel- 
atively bright (J<16) CH^BDs in 1800^ sq. degrees and 
the three much fainter candidates discussed here. 

We also showed that, due to their very low intrinsic 
luminosities, methane dwarfs become numerous (relative 
to L and M dwarfs) only at very faint J magnitudes (see 
figures from |^ to |6|). To significantly improve the ob- 
servational database it is therefore necessary to cover a 
reasonably large area of sky (several sq. degrees) down to 
J~22, as well as K~22 and 1^27, the latter two being 
necessary for selecting the CH4BD candidates. We look 
forward for this type of observations to be soon performed 
by the IR large field deep imagers operating or planned 
on large telescopes. 
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